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ABSTRACT 

Using a sample of 28 H II regions from the literature with measured tem- 
perature inhomogeneity parameter, t"^, we present a statistical correction to the 
chemical abundances determined with the Te (4363/5007) method. We used the t"^ 
values to correct the oxygen gaseous abundances and consider the oxygen deple- 
tion into dust to calculate the total abundances for these objects. This correction 
is used to obtain a new calibration of Pagel's strong-line method, R23, to deter- 
mine oxygen abundances in H II regions. Our new calibration simultaneously 
considers the temperature structure, the ionization structure, and the fraction 
of oxygen depleted into dust grains. Previous calibrations in the literature have 
included one or two of these factors; this is the first time all three are taken 
into account. This recalibration conciliates the systematic differences among the 
temperatures found from different methods. We find that the total correction 
due to thermal inhomogeneities and dust depletion amounts to an increase in the 
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0/H ratio of H II regions by factors of 1.7 to 2.2 (or 0.22 to 0.35 dex). This 
result has important imphcations in various areas of astrophysics such as the 
study of the higher end of the initial mass function, the star formation rate, and 
the mass-metallicity relation of galaxies, among others. 

Subject headings: galaxies: abundances - H II regions galaxies: ISM - H II 
regions regions - H II regions ISM: abundances 



Introduction 



The study of physical conditions of H II regions is a fundamental tool to constrain 
models of nucleosynthesis of massive stars, galactic chemical evolution, and the chemical 
evolution of the Universe. In the same way, the study of the chemical composition of metal- 
poor H II regions is crucial for the understanding of the primordial chemical composition 
of the Universe. Most of these studies have been done using abundances obtained from the 
so called direct method, which assumes a homogeneous temperature throughout the whole 
volume of the object and then uses this temperature to determine abundances of all available 
ions; this temperature is commonly obtained from the [O III] coUisionally excited line ratio 
AA 4363/(4959+5007); hereafter we will call this procedure the Te(4363/5007) method. 

Unfortunately, the simplification of the thermal structure of H II regions implied by the 
Tg (4363/5007) method, yields lower abundances than those obtained with other methods. 
For example: (i) the 0/H abundances derived from recombination lines, (ii) the 0/H abun- 
dances derived by comparing the nebular lines with photoionization models, and (iii) the 
comparison of the 0/H nebular abundanc es (including gas and dust), with those derived 
from B stars and young F and G stars (e.g. 



Przvbilla et al.ll2008l : iPeimbert fc Peimbert 



Bensbv fc Feltzindl2006l: 



2011 



Kew. 



Simon-Dfaz fc Stasihska 



ev fc Elhson 20081: 



2011 



iNieva. fc Przybilla 



2OI2I ). The abundances derived from these three different methods are consistent within the 
errors. The abundances obtained through coUisionally excited lines corrected for the presence 
of thermal inhomogeneities are consi stent with the abundan c es ob taine d through other meth 



ods. For reviews on such issue see iPeimbert fc PeimbertI (I2OIII ) and iLopez-Sanchez et al. 
J2012h . 



Abundances obtained with the Te (4363/5007) method are underestimated because the 
temperature used in the determination process is assumed to be uniform. In the presence of 
either local or global thermal inhomogeneities, coUisionally excited lines will be brighter in 
the hot regions while recombination lines will be brighter in the cool regions. This leads to 
choosing a temperature higher than the average and thus to assume that less oxygen (and 
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more hydrogen), than the amount that is actually present, is needed in order to reproduce 
the observed line intensities. For this reason, these abundances should be corrected for the 
presence of thermal inhomogeneities before being used to calibrate any strong-line indicator. 



Pagel et al.l ( 1l979l ) were the firsts to introduce the i?23 versus 12-|-log(0/H) diagram 
(from here on, this diagram will be referred to as the i?23^0/H diagram). They found 
that the sum of the intensities of nebular oxygen with respect to H/3, J([0 II] 3727) + 
/([O III] 4959+5007) /J(H/3), varied smoothly with the total oxygen abundance and pro- 
posed to use what is now known as i?23, O23, or Pagel's method to determine the chemical 
abundances of objects where no auroral lines could be observed and, hence, no tempera- 
ture could be determined. At present, more than a dozen strong- li ne indicators have been 



proposed; of them, the rn ost studied is R22, (e.g. iStasihska] l2006l : iKewley fc EUisonl 12008 
Lopez-Sanchez et al.ll2012l ). These strong-line indicators can be calibrated using well known 
abundance determinations such as the ones based on the Te (4363/5007) method or from 
fitting lines with numerical photoionization models. 

The quahty of the abundances derived with the Te (4363/5007) method, strongly de- 
pends on the validity of the assumption of constant temperature over the observed value. In 
the presence of large thermal inhomogeneities, the Tg (4363/5007) method will underestimate 
abundances since the temperature dependence of the intensity of collisionally excited lines 
(such as [O III] AA 4363, 4959, and 5007) is very different from that of the recombination line 
H/9. The intensity of collisionally excited lines is proportional to / oc Te"^''^ xexp(— A E/YT^), 
while the intensity of H/3 is proportional to / a Tg~°'^'', where Te is the local electron temper- 
ature, k is the Boltzmann constant, and A E is the collisional energy required to excite the 
line. As an alternative to avoid the problem of the large temperature dependence of collision- 
ally excited lines, recombination lines can be used to determine abundances and calibrate 
the i?23^C)/H diagram because they have a temperature dependence given approximately 
by / oc T^^. For recombination line to H/9 ratios, the temperature dependence almost 
cancels out, whereas for collisionally excited line to H/3 ratios the temperature dependence 
does not cancel out. Some work has been done in the upper branch of the diagram using 
recombination lines, however, these faint lines can only be observed in metal-rich objects. 

The discrepancy between abundances determined using recombination lines and colli- 
sionally excited lines is called the abundan ce discrepancy factor , ADF^ problem. Typical val- 
ues for H II regions lie in the 1.5 to 3 range ( 



2007; Peimbert et al. 2007; Esteban et al. 



T ^eimbert et al.lll993l . l2005l : iGarcia-Rojas fc Esteban 
20091 and references therein). There have been 



two major explanations for this discrepancy: (i) high metallicity inclu sions that will create 



cool high density regions surrounded by hot low density regions (e.g. iTsamis fc Pequignot 



20051 ) , and (ii) thermal inhomogeneities in a chemically homogeneous medium that are caused 
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by various physical processes such as sha dowed regions, advancing io nization fronts, shock 



waves, magnetic reconnection, etc. (e.g. iPeimbert &: PeimbertI l201ll ). Since the densities 



determined from oxygen rec ombination hues correlate very well with the [CI III] densities 
( IPeimbert fc PeimbertI |2005[ ). we prefer the second explanation. 



In Section |2] we discuss the corrections due to thermal inhomogeneities and due to the 
fraction of O depleted into dust grains. In Section |3] we present a quantitative correction 
to the Tg (4363/5007) method, this correction is applied to the i?23 niethod in Section HI In 
Section O we discuss the relation between the 0+"'" fraction and the intensity of the [O II] 
and [O III] nebular lines. Finally, in Section [6] we present the conclusions. 



2. Temperature Inhomogeneities and Dust 



The formalism to study the chem ic al ab undan ces in the presence of therm al inhomo- 
geneities was introduced by iPeimbertI ( 1l967l ) and IPeimbert fc Costerd (119691 ) . where the 
parameter that characterizes such thermal inhomogeneities was defined as t^. In these works 
it was recognized that thermal inhomogeneities would produce systematic differences in tem- 
peratures determined from different methods, and that abundances determined with each 
one of those temperatures would also present systematic differences. The abundances deter- 
mined using the Te(4363/5007) method were approximately a factor of 2.5 lower than those 
determined using the same coUsionally excited line intensities including the effect of thermal 
inhomogeneities. 



There are several ways to obtain a value, six of them are discussed in lPeimbert fc Peimbert 

( I2OIII ). These are based on the comparison of: (i) the O"*""*" abundances obtained using [O III] 



lines versus abundances obtained using O II lines; (ii) the temperature derived from He I 
lines versus the temperatures derived from [O II] and [O III] lines; (iii) the temperatures 
derived from the ratio of the Balmer and Paschen continua to the Balmer line intensities 
versus the temperatures derived from [O II] and [O III] lines; (iv) the C"*"^ abundances de- 
rived using [C III] and C III] lines versus abundances obtained from C II lines; (v) 1.5x10^ 
columnar teraperat ures determined from a high spatial resolution map of the Orion nebula 
( lO'Dell et al.ll2003l ): and (vi) the calibrations of the i?23 method using observations versus 
those using photoionization models. In all cases, the values are consistent with each other, 
which implies that those objects are chemically homogeneous. The first two methods are the 
most used to determine t^. 



Independently of the effects of the thermal structure, there is an additional correc- 
tion due to dust depletion that should be taken into account to determine abundances. 
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Esteban et al.l (119981 ) found that the depletion of oxygen into dust grains in the Orion nebula 



amounted to a 0.08 dex correction in the 0/H ratio; slightly higher estimates of this cor- 



rection have been found in re cent studies of the Orion nebula (IMesa-Delgado et al.l l2009b 



Simon-Diaz fc Stasihskal l201l[). From the dep l etion of Mg, Si, and Fe in Galactic and ex- 
tragalactic H II regions, IPeimbert fc PeimbertI J2010h estimated an 0/H correction ranging 
from 0.09 to 0.11 dex that increases with increasing 0/H ratios. We have used the correction 
recommended by Peimbert et al. for all the objects of our sample; this corrections amounts to 
0.11 dex for those objects whose gaseous oxygen abundance is 12+log(0/H)> 8.3, 0.10 dex 
for objects with 8.3 >12+log(0/H)> 7.8, and 0.09 for those objects with 7.8 >12+log(0/H). 



3. Correction to the abundances derived with the Te(4363/5007) method 



Due to the inconsistencies presented when comparing abundances obtained with the 
Te (4363/5007) method and those obtained with, recombination lines, photoionization mod- 
els, and using abundances of recently formed stars; it is clear that abundances derived from 
the Te(4363/5007) method are systematically underestimated. This needs to be corrected 
by considering the thermal structure of the object and the fraction of oxygen depleted into 
dust grains. 

Table [T] lists the sample we used to correct the Tg (4363/5007) method. The sample 
consists on six H II galaxies, eight Galactic H II regions, and 14 extragalactic H II regions. 
This total of 28 objects were gathered from the literature and, in each of them, at least one 
value of was obtained. For the cases where there was more than one value we took the 
one presented as average or the one from the brightest observed region. 

Based on detailed analysis presented in the literature it follows that the total oxy- 
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|2]). We expect the magnitude of these effects to depend on the oxygen abundance and, 
possibly, on the oxygen excitation degree. In order to take these factors into account, we 
considered that the abundances derived from the Tg (4363/5007) method should be corrected 
with a function / dependent on the oxygen ionization degree and on the oxygen abundance 
determined from the Te(4363/5007) method, (0/H) 4363/5007: 



(0/H)cALM = (O/H)4363/5007 + /(-P, (O/H)4363/5007) , 



(1) 



- 6 - 



where CALM stands for Corrected Auroral Line Method. Throughout this paper all oxygen 
abundances are given in units of 12+log(0/H). 

As a first approximation to the function /, we decided to use a linear correction of the 

form 

/(P, (O/H)4363/5007) = + C2 X (O/H)4363/5007 + X P + X P X (O/H)4363/5007, (2) 

where Ci, C2, C3, and C4 are constants. 

All abundances determined for our sample were corrected for (i) the presence of ther- 
mal inhomogeneities and (ii) the fraction of oxygen depleted into dust grains. We have 
found that the difference between these corrected abundances and those determined with 
the Te(4363/5007) method lie in the 0.16 to 0.40 dex range. 

We found very little dependence with oxygen degree of ionization: C3 and C4 are 0.010 
and 0.0002, respectively. Since the dependence with ionization degree is small and consistent 
with zero we will ignore it. Figure [1] presents the oxygen abundance correction dependence 
with oxygen abundance for the objects of our sample. The best fit is obtained with Ci and 
C2 equal to —0.375 and 0.0825, respectively: 

/(P, (O/H)4363/5007) = /((O/H)4363/5007) = 0.08 25(O/H)4363/5007 " 0.375; (3) 

this implies a corrected oxygen abundance: 

(0/H)cALM = 1.0825(O/H)4363/5007 " 0.375. (4) 

The points in Figure [1] show considerable dispersion; this is not surprising, because (i) the 
errors in the determination of each t"^ value are large and (i i) the true value of for each 



object depends on the details of the object (see Figure 4 of iPeimbert et al.l 120121 ). To test 
for the statistical significance of this correlation we calculated a Spearman's correlation 
statistic of p = 0.3685, this suggests that the correlation is real with 95% confidence. When 
comparing the dispersion with the uncertainties on the oxygen abundance determinations 
for this sample we calculated = 60.28; showing that about half of this dispersion comes 
from the uncertainties of the determination of t^; meaning that approximately half of the 
dispersion comes from the unique physical conditions of each object. 



4. Recalibration of the P23 Abundance Determination Method 



R23 is now being widely used to determine abundances in objects with low intrinsic 
brightness or objects with low to intermediate redshift (2; up to 1). The P23 method remains 
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one of the most popular strong line indicators because (i) it is based directly on oxygen lines, 
(ii) it includes lines from the 2 relevant ionization stages, (iii) the lines it uses are the most 
intense for most H II regions, and (iv) these lines lie in the blue part of the spectrum and 
can be observed from the ground for higher redshifts than other methods. 

Unfortunately, the curve in the i?23^0/H diagram is double valued (see Figure [2]); this 
occurs because: at low metallicities, the main cooler in H II regions is the hydrogen atom; 
with increasing element abundances, the cooling due to oxygen nebular lines increases, hence, 
the value of R23 increases reaching a maximum value of R23 ~ 10 when 12+log(0/H)isM ~ 
8.4; then, at higher metallicities, the IR [O III], IR [Ne II], and optical [N II] lines begin 
to dominate the cool ing so that the total value of R23 decreases as the metallicity increases 
dVan Zee et"aDll998l ). 



Also, there is a noticeab le dependence on the oxygen ionization de gree fe.g. iMcGaugh 



1991 



Kewley fc Dopitall2002[ l. In a series of articles iPilyugiru (120001120011 ) and lPilyugin &: ThuanI 



fcoOS ). presented an analytic family of curves for each branch of the R23 - 0/H diagram. This 
family of curves depends on the oxygen excitation degree across the object; defining the P 
parameter as the observable oxygen excitation ratio, P=([0 III] 4959+5007)/([O II] 3727 + 
[O III] 4959+5007), and presenting two equations to determine the upper and lower branch 
abundances as a function of (0/H)p=/(P,i?23)- This calibration is based on abundances 
determined using the Tg (4363/5007) method. 

Unfor tunately, depending on the calibration, result ant O /H values can vary by as much 
as 0.6 dex ( Peimbert et al.|[2007 : Kewley fc Ellison 2008). This large dispersion in the calibra- 
tions mainly comes from what was chosen to be adjusted: (i) the theoretical intensities given 
by photoionization models, (ii) abundances derived from the Te (4363/5007) method, or (iii) 
abundances derived from observations of oxygen recombination lines. The lowest determina- 
tions for the gaseous oxygen abundances are those that come from using the Te (4363/5007) 
method, therefore they require the correction due to thermal inhomogeneities; additionally, 
to obtain the total oxygen abundances, a correction of ~ 0.10 dex — due O depletion into 
dust grains — is required. 

We find that, among all the calibrations of the R 23-O/B. diagram available in the 
literature, the one presented by lPilyugin fc ThuanI (120051 ) is the most adequate for this work 
since it explicitly includes the oxygen ionization degree through the P parameter. 



The calibration of iPilyugin fc ThuanI ( l2005l ). (0/H)p, used abundances that were de- 
termined with the Te (4363/5007) method; hence, abundances are underestimated but this 
allows us to use equationlHto correct them since (0/H)p has abundances determined with the 
Te (4363/5007) method at its core. We can determine oxygen abundances from a recalibration 
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of the i?23 method, (0/H)rrm, by substituting (0/H)rrm and (0/H)p in place of (0/H)calm 
and (0/H)43e3/5oo7 into equation (jl]). We take (0/H)p from equation [22] and [24] of Pilyugin 
and Thuan for the upper and lower branches, respectively, to obtain two recalibrated families 
of curves that now simultaneously consider the thermal structure, the ionization structure, 
and the depleted fraction of O into dust grains. The oxygen abundance for the upper branch 
(which shifted from 12+log(0/H)p >8.25 up to 12+log(0/H)RRM >8.55), can be estimated 
from: 



(0/H)™m = I fi^^. + 7261 + 8 42P + 327,5P^ V.0825 - 0.375 

V / ^HHM \85.96 + 82.76P + 43.98P2 + 1.793/^23 y 
i?23 + 1837 + 2146P + 850P2 

(5) 



209.5 + 201. 7P + 107.2P2 + 4.37P. 



23 



while the oxygen abundance for the lower branch (which shifted from 12+log(0/H)p <8.00 
up to 12+log(0/H)RRM <8.29), can be estimated from: 



, , R23 + 106.4 + 106.8P - 3.40P2 \ 

'O/"'™" = ' l7.72 + 6,60P + 6.95P^-0,302flJ ^°«^^-°'^^^ 
P23 + 90.73 + 94.58P - 5.26P2 

(6) 



14.81 + 5.52P + 5.81P2 - 0.252^23 



Figure |2] shows these recalibrations of the upper and lower branches for different oxygen 
io nization degrees, w h ile Fi gure [3] compares the calibration presented in this work and that 
of iPilyugin fc ThuanI ( 120051 ) to the data from Table [TJ 



5. The 0++ fraction 

Something missing when using the P parameter is a quantitative connection with the 
O^"^ fraction in the nebulae. We will define the O^"*" ionization degree as: 

OID = /^en(0++)cir 

J nMO+ + 0++)dr ^ ' 

Since there is a tight relationship between the parameter P and the O^'^ ionization 
degree (see Figure H]), the fraction of 0+"'' in the object can be estimated through the the 
following equation: 

OID = 0.9821P - 0.0048. (8) 
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This equation allows us to statistically estimate the O"'""'" ionization degree directly from the 
observable P without having to make a lengthy analysis of the ionic abundances. It should 
be noted that while this relation shows no significant bias with the 0/H ratio, more points 
are needed to confirm the lack of bias or to estimate the relevance of the oxygen abundance 
in this correlation. 



6. Conclusions 

We used a sample of 28 H II regions obtained from the literature, with values ranging 
from 0.020 to 0.120, to present an important correction to the abundances derived with 
the Te(4363/5007) method, as well as a recalibration of Pagel's strong-line method, -R23 to 
determine chemical abundances. Both, the correction and the recalibration consider the 
thermal and ionization structure of the object, as well as the fraction of oxygen depleted 
into dust grains. We found that the differences between the corrected abundances and those 
determined with the Te(4363/5007) method lie in the 0.22 to 0.35 dex range. 

Almost all previous calibrations of the -R23^0/H diagram in the literature have con- 
sidered a homogeneous temperature throughout the object. The recalibration presented in 
equations O and O is the first to simultaneously consider the thermal structure of the pho- 
toionized object, the ionization structure, and the correction due to the fraction of oxygen 
depleted into dust grains. This recalibration is particularly useful for objects with low to 
intermediate redshifts (up to z ~ 1). 

We find that the O^"*" ionization degree shows a surprisingly almost one-to-one rela- 
tion with the parameter P, thus ensuring that a calibration which considers P is directly 
considering the oxygen ionization structure of the object. 

Given that all strong-line indicators attempt to determine the total oxygen abundance 
as means of estimating the heavy element content, and that all the observational calibrations 
have been done using abundances determined with the Te (4363/5007) method, the oxygen 
correction presented in this paper should be used to recalibrate other strong-line indicators 
in the literature. 

The up-shifting of the family of curves in the -R23^0/H diagram found in this work, 
has important consequences in various areas of astrophysics such as the determination of: 
initial mass functions, star formation rates, luminosity functions, yields, radial metallicity 
gradients of spiral galaxies, and chemical evolution models of the Universe. 

We are grateful to an anonymous referee for a careful reading of the manuscript and 
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Table 1. H II region data 
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NGC 3576 


GR 


8500±50 


0.038±0.009 


0.78 


0.78 


0.67 


8.56 


8.92 


1 


M16 


GR 


7650±250 


0.039±0.006 


0.59 


0.28 


0.25 


8.50 


8.90 


2 


M17 


GR 


8950±380 


0.033±0.005 


0.73 


0.83 


0.83 


8.52 


8.88 


3 


M8 


GR 


8090±140 


0.040±0.004 


0.53 


0.38 


0.28 


8.51 


8.85 


3 


H1013 


XR 


7370±630 


0.037 


0.42 


0.49 


0.49 


8.45 


8.84 


4 


NGC 595 


XR 


7450±330 


0.036 


0.51 


0.37 


0.44 


8.45 


8.80 


4 


M20 


GR 


7800±300 


0.029±0.007 


0.60 


0.20 


0.17 


8.53 


8.79 


2 


Orion 


GR 


8300±40 


0.028±0.006 


0.77 


0.86 


0.83 


8.51 


8.79 


5, 6 


NGC 3603 


GR 


9060±200 


0.040±0.008 


0.89 


0.92 


0.93 


8.46 


8.78 


2 


K932 


XR 


8360±150 


0.033 


0.72 


0.72 


0.79 


8.41 


8.73 


4 


NGC 2403 


XR 


8270±210 


0.039 


0.59 


0.66 


0.67 


8.36 


8.72 


4 


NGC 604 


XR 


8150±160 


0.034±0.015 


0.60 


0.71 


0.71 


8.38 


8.71 


4 


S 311 


GR 


9000±200 


0.038±0.007 


0.72 


0.32 


0.31 


8.39 


8.67 


7 


NGC 5447 


XR 


9280±180 


0.032 


0.85 


0.78 


0.86 


8.35 


8.63 


4 


30 Doradus 


XR 


9950±60 


0.033±0.005 


0.90 


0.85 


0.85 


8.33 


8.61 


8 


NGC 5461 


XR 


8470±200 


0.027±0.01z 


0.71 


0.80 


0.77 


8.41 


8.60 


4, 9 


NGC 5253 


HUG 


11960±290 


0.072±0.027 


0.96 


0.85 


0.78 


8.18 


8.56 


10 


Transition zone 




















NGC 6822-V 


XR 


13000±1000 


0.076±0.018 


0.91 


0.88 


0.89 


8.08 


8.45 


11 


NGC 5471 


XR 


14100±300 


0.082±0.030 


0.93 


0.75 


0.78 


8.03 


8.33 


9 


NGC 456 


XR 


12165±200 


0.067±0.013 


0.83 


0.78 


0.80 


7.99 


8.33 


12 


Lower branch 




















NGC 346 


XR 


13070±50 


0.022±0.008 


0.92 


0.88 


0.69 


8.07 


8.23 


13, 14 



Table 1 — Continued 



Object 


Type'^ 


Te[0 III] 




log(i?23) 


P 


0++ 
0++0++ 


O/H^^ 


0/H^ 


References 


NGC 460 


XR 


12400±450 


0.032±0.032 


0.81 


0.62 


0.56 


7.96 


8.19 


12 


NGC 2363 


XR 


16200±300 


0.120±0.010 


0.92 


0.97 


0.97 


7.76 


8.14 


4 


TOL 2146 - 391 


HUG 


15800±170 


0.107±0.034 


0.91 


0.92 


0.86 


7.79 


8.09 


15 


TOL 0357 - 9315 


HUG 


14870±230 


0.029±0.064 


0.93 


0.93 


0.87 


7.90 


8.12 


15 


Haro 29 


HUG 


16050±100 


0.019±0.007 


0.91 


0.91 


0.88 


7.87 


8.05 


13, 16 


SBS 0335-052 


HUG 


20500±200 


0.021±0.007 


0.67 


0.93 


0.93 


7.35 


7.60 


13, 17 


I Zw 18 


HUG 


19060±610 


0.024±0.006 


0.47 


0.86 


0.90 


7.22 


7.41 


13, 17 



(1200911: f5 ^ 



Peimbert et al 



Esteban et al 



2004|lli6l 



Refere nces.— fl ) lGarcfa-R,o1as et all (12004 h ('2') iGarcia-Eoias et al l mm-. ('3'! iGarcia-Roias et al l mm-. 



O'Dell et al 



(l2ooii ■ f7HGarcia-Roias et al 



Esteban et al] (l2002l ): flO ) iLopez-Sanchez et al.l (|2007l); (11) IPeimbert et al.l (|2005I M12} 



200. 



Esteban et al. 
PeimbertJ (120031'): ( 9'! 



Pena- Guerrero et al 



i2012l): fl3) 



\2mi \. (14) IPeimbert et al.l (|20od ): (15) IPeimbert etaP (|2012l ): (16) llzotov et al.l (Il997l ): (17) llzotov fc Thuan 

(|l999l ). 



'^GR=Galactic region, XR=Extragalactic region, HIIG= H II Galaxy. 

Gaseous O abundance with homogeneous temperature, 4^=0.000. In units of 12+log(0/H). 

^Total O abundance with thermal inhomogeneities, t^>0.000, including the correction due to depletion of O into 
dust grains. In units of 12+log(0/II). 
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7 7.5 8 8.5 

12 + lOg(0/H)4363/5007 



Fig. 1. — Plot of the abundances obtained with the Tg (4363/5007) method, versus 
the difference in dex between the corrected abundances and those obtained with the 
Te(4363/5007) method. The hne is the best hnear fit to the objects presented in Table 

m log(O/H)iSM-log(O/H4363/5007)= 0.08 25x[12+log(O/H4363/5007)]-0.375. 
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Fig. 2. — New calibration of the R23 Method from equations [5] and [61 For clarity, the solid 
P-valued curves have been labeled. The dashed horizontal lines indicate the transition zone 
between the upper and lower branches. 
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Fig. 3. — Plot of the oxygen abundance vs. -R23 for the 28 objects from Table [H where 
abundances were determined considering the presence of thermal inhomogeneities as well as 
the fraction of O depleted into dust grains. For comparison we present representative bands 
of the upper and lower branches for both nebular line calibrations, (0/H)rrm (solid black 
curves) and (0/H)p (point-dashed curves). 
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Fig. 4. — Linear relation found between the observational parameter P and the true O 
ionization degree. The objects used to obtain this relation are those listed in Table [H 



